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\mJ • The account of the origin of cosmic structure, as provided by the standard inflationary paradigm, 

Q^i is not fully satisfactory, as has been argued in [A. Perez, H. Sahlmann, and D. Sudarsky, Class. 

2 ' Quantum Grav., 23, 2317, (2006)]. The central point of that work is to point out the need to discuss 

^ ' and explore the physical mechanism that is capable of generating the inhomogeneity and anisotropy 

of our Universe, starting from an exactly homogeneous and isotropic initial state associated with 
^ I the early inflationary regime. We review this issue briefly here together with the proposal to 
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address this shortcoming in terms of a dynamical collapse of the vacuum state of the inflaton field. 
We also briefly indicate how this issues might be connected to other questions being faced in the 
^^ , study of the quantum/gravity interface, and their relevance to the investigations concerning the 



statistical characterization of the primordial spectrum. 
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I. INTRODUCTION 

At what point in the cosmic evolution do the actual primordial inhomogeneities arise?. I.e 
when does our universe depart from the homogeneity and isotropy that results from inflation. 
This is a question that one might expect should be addressed, at least in principle, by any 
theory which deals with the emergence of cosmic structure. Yet, in the standard inflationary 
account [Ij, which is nowadays regarded as a remarkable success, the context where such 
issues can be addressed seems to be simply absent pi]. That is, within the orthodox accounts, 
one can not identify the physical process responsible for the generation of those features in 
our Universe. In fact, according to the inflationary paradigm, from a relatively wide initial 
set of possibilities marking the end of the mysterious quantum gravity era, the accelerated 
inflationary burst leads to a homogenous and isotropic (H&I) Universe where the quantum 
fields are all characterized by the equally homogeneous and isotropic vacuum states (usually 
taken specifically to be the so called Bunch-Davies vacuum). From these conditions, it 
is usually argued, in a rather uncleaiu although strongly image-evoking manner, that the 
quantum fluctuations present in such quantum state morph into the seeds of anisotropics and 
inhomogeneities that characterize our late Universe. This issue is sometimes characterized 
as the "transition from the quantum regime to the classical regime" , but we flnd this a bit 
misleading: most people would agree that there exist no distinct and separated classical and 
quantum regimes. The fundamental description ought to be always quantum mechanical; the 
so called classical regimes are those where certain quantities can be described to a sufficient 
accuracy by their classical counterparts representing the corresponding quantum expectation 
values. The paradigmatic example here is provided by the coherent states of an harmonic 
oscillator which correspond to minimal wave-packets with expectation values of position 
and momentum following the classical equations of motion. In any case, it seems clear that 
from a situation corresponding to a H&I background and a quantum aspect characterized 
by a H&I state, one can not end upcl in a situation which is characterized, at the classical 
level, as containing actual inhomogeneities and anisotropics. It is clear that, in terms of 
the dynamics, such transition can not be accounted for by the gravity/inflaton action which 



^ Acknowledgments that this is an unclear aspect can be seen for instance in Weinberg's Cosmology 

where the author explicitly states his view on the subject. 
^ In the absence of something else, which in other circumstances would be identified as a measurement, but 

which clearly can not be invoked in the present setting. Observers and measuring apparatuses are only 

possible well after the H&I has been broken, so t^se can hardly be part of the cause of the breakdown. 



preserves such symmetries. Simply put, if the initial state is H&I and the Schrodinger 

evolution is tied to a Hamiltonian that preserves these symmetries, the resulting state can 

not be anything but H&I. Nonetheless, various types of arguments are often put forward 

in attempts to bypass the above conclusion. The discussion of the conceptual problems 

and those associated have been discussed in previous works by some of us and by others in 

4-|6|. We will not reproduce those arguments here as the objective of the present work is to 

focus in the statistical aspects that emerge tied to what we consider as a more conceptually 

transparent picture, and which is based on a modified version of the standard inflationary 

paradigm, that we have been advocating in previous works [4, [Zl, |8| . 

We should note however, that we can not escape from the related problems even if we 

choose to adopt a very pragmatic position. Lets assume that one chooses to ignore the 

shortcomings of the more standard accounts, and accepts that, say decoherence, addresses 

somehow, the issue at hand, and that the mystery lies only in the question concerning the 

precise mechanism that lies behind the fact that, from the options exhibited in those analyses 

(i.e. those displayed in the reduced density matrix), one single particular option seems to 

be selected J9| . Within this point of view, one would be assuming that the initial symmetry 

has been lost -at least for practical purposes, as presumably one would be advocating 

when adopting such position- in association with that particular realization or actualization 

(represented by a particular element in the density matrix). Thus, it seems clear that for the 

sake of self consistency, one should consider, when studying aspects of the inhomogeneity 

and anisotropics in the CMB, the state corresponding to such "selected option" and not the 

entire vacuum state which describes the H&I state of affairs previous to the "selection" u. 

In following such views, the discussion that we are presenting in this paper, would have to 

be taken to represent the effective description corresponding to "our perceived Universe" (in 

a context where one puts together something like the many-worlds interpretation, with the 

arguments based on decoherence). Although we definitely do not adhere such view for the 

reasons explained in [5[, it is clear that an effective description such as the one presented 

here is what would have to be contemplated when dealing with the issues within any view 

^ The selection of course refers to the fact that according to the standard arguments the resulting density 

matrix, after becoming essentially diagonal due to decoherence, represents an ensemble of universes and 

our particular one corresponds to one of them. That one can be considered as selected by nature to 

become realized. Alternatively, one might take the view that these other universes are also realized, and 

thus they also exits in realms completely inaccessible to us. In that case the selection corresponds to that 
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which pretends to allow one deal with the details characterizing the inhomogeneities and 
anisotropics in the cosmic structure and its imprints in the CMB that we do observe. 

In fact, it is perhaps worth noting, that in the usual accounts, it is hard to pinpoint 
where exactly the statistical aspects come into play and which kind of the statistics one 
is dealing with. That is, in the standard approach, the specific Universe is not described 
in any sense (not even in terms of any unknown, yet explicitly identified quantities), and 
the randomness which characterizes it lies hidden in unspecified aspects associated with the 
multiple interpretations. In other words, one can not identify the random variables; one does 
not know how many there are; one can not say how the various elements of the ensemble 
differ from each other, etc. 

In order to fully and satisfactorily address the problem, it seems paramount to be able 
to point out what exactly is wrong with the argument leading to the conclusion of the the- 
orem drawn above, i.e. where does nature deviate from the theory leading to the erroneous 
conclusion that our Universe is even today, at the fundamental quantum level, perfectly ho- 
mogeneous and isotropic?. It follows that such explanation must indicate where, the ordinary 
^/-evolution -with the symmetry preserving Hamiltonian- does break down. We can easily 
see that none of the proposals to deal with the issue, based on the standard paradigms, can 
point to any reason where that breakdown might occur, and much less, point to a physical 
reason for that departure from standard Quantum Theory. This has led us to take a view to 
tie this problem with the ideas advocated by Roger Penrose, who argues [10[ that quantum 
theory should itself suffer modifications as a result of its combination with the fundamental 
theory of space-time structurcl. Among the aspects of the theory that would be substan- 
tially affected according to Penrose's views, are those related to the reduction postulate (or 
R process, as he calls it), and its contraposition -and to some degree its contradiction- with 
the unitary evolution (or U process, as labeled in his works) controlled by Schrodinger's 
equation. In fact, the issue of dynamical quantum reduction has received a lot of attention 
within the community working in foundational aspects of quantum theory, and there are in 
the existing literature several well defined proposals such as those in lll-ll3|. 

The proposal behind our work is based on the inclusion of the hypothesis that a dynamical 

'^ This is what is often though as quantum gravity. We did not use that term because that often presupposes 

that one is considering the relevant theory to be simply the adaptation of general relativity to the standard 

quantum theory, while what one has in mind when following Penrose's ideas is something much more 

distant from known physics, involving as indicated, modifications of quantum theory itself. 
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collapse of the wave function lies behind the breakdown of the initial homogeneity and 
isotropy. In other words, that a non- unitary jump in the quantum state plays a role in 
transforming the inflaton vacuum into a quantum state that lacks the translational and 
rotational symmetries of the former state. 

It goes without saying that we can not at this stage try or hope to point out the physical 
origin of such dynamical collapse. However, once one has accepted that something of this sort 
is occurring, one can parametrize its basic characteristics, and use the relevant observational 
data to infer something about the nature of the new physics that lies behind the phenomena. 
This has been the basic attitude behind the program stared in [^. 

In this paper our aim is to briefly explore for the first time some of the basic differences 
associated with the statistical considerations, between the usual account and the proposal 
we have been working with. 

In order to be a bit more explicit lets start by reminding the reader that in the standard 
approaches, the study of the statistical nature of the problem is based on the study of the 



ra-point functions of the Newtonian potential, \l/(xi) . . .\l/(x„), with the overline denoting 
the average over an ensemble of Universes. There one needs to face the issue of what is the 
relationship between the quantum n-point functions and the quantities we actually measure. 
The usual approach is based on the identification of 



(0|^(xi) . . . *(a;„)|0) = $(xi) . . . ^(a;,), (1) 

that is, based on the identification of quantum and statistical n-point functions. The latter 
are naturally associated with ensembles of Universes, all of which, even if real, are unacces- 
sible to us. Therefore, invoking ergodic arguments, a further connection is made between 
ensemble averages and time -and in our case- spatial averages. 

At this point we must express some caution regarding such identifications and to critically 
explore the validity of many of underlying assumptions, and also to unearth the places where 
the issues at hand can have important statistical/observational consequences. 

In order to further motivate this work it is perhaps worth pointing out the simplest 
place where we can appreciate the problematic aspects of unquestioningly accepting such 
identifications in a very simple situation concerning the one point function: Consider the 



standard treatments based in the so called Mukhanov-Sasaki variable: 

where \E' is the metric perturbation known as the Newtonian potential, 0o is the derivative 
of the background inflaton with respect to conformal time rj, 6(j) is the perturbation in the 
inflaton field, a is the scale factor and "H = - (related to the standard Hubble parameter H 
through "H = aH). The Einstein equations then lead to Au = -2 (f ) and v = ^ i^u)' where 
z = ^. Given the equations of motion, the Newtonian potential can thus be expressed in 
terms of the field v{x,ri) and its momentum conjugate 'iT^{x,ri) = v{x,ri). The expression 
for the corresponding Fourier components is 

nM = -^^ {-,,(,) - ^M,)) . (3) 

where we e is the so-called slow-roll parameter e = 1 — 'H/'H?. 

We are interested in the temperature anisotropies of the CMB observed today on the 
celestial two-sphere, which are related to the inhomogeneities in the Newtonian potential on 
the last scattering surface, 

^(^,V^) = ^^(^D,fD). (4) 

The data is described in terms of the coefficients aim of the multipolar series expansion 
^(^,V^) = 5^«,„r,„(^,V^), aim= [^^iO,^)Y^liO,^)d^- (5) 

Jo , J -Lq 

Im 

Here 6 and ip are the coordinates on the celestial two-sphere, with F/m(^, V^) the spherical 
harmonics. 

The quantities aim are then given by 

«/™ = ^ y j^s3i{kRD)Y,l{k)A{k)^^,M, (6) 

with jilkRc) the spherical Bessel function of order /; rjji is the conformal time of reheating 
which can be associated with the end of the infiationary regime and Rd the comoving radius 
of the last-scattering surface. We have explicitly included the modifications associated with 
late time physics encoded in the transfer functions A{k). 

Now, the problem is that if we compute the expectation value of the left hand side (i.e. 
identifying (\l/) = \1') in the vacuum state |0) we obtain 0, while it is clear that for any given 
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/, m the measured value of this quantity is not 0. That is, if we rely in this case on the one- 
point function and the standard identification, we find a large conflict between expectation 
and observation. We might even be tempted to say that evidence of non-Gaussianity has 
already been observed in each measurement of a particular aim- This is, of course, not 
what one wants, however this makes clear that disentangling the various statistical aspects 
(ensemble statistics; space and time statistics, including orientation statistics and finally 
the nature of the assumed connection of quantum and statistical aspects) and to make 
explicit the assumptions underlying the identifications, as well as the expected limitations, 
is paramount to avoid the dangers of confusion. 

As a matter of fact, what exactly is wrong with the above argument?. Normally one 
would have to argue that the relevant prediction concerns the ensemble averages and thus 
identifying the result 

^im = ^J ^,Ji{kRD)YCm{k)A{km^^,M\0) = 0, (7) 

with the estimate of the actually measured quantity aim would be incorrect. One could 
however wonder why would it be incorrect to say that we might rely on ergodic arguments 
and thus expect that spatial averages would correspond to the ensemble averages, and note 
that the measured quantity corresponding to the expression above is after all already a 
weighted average over the celestial two sphere. The answer would need to be that we 
still must carry out a further average: the average over orientations in order to have any 
confidence that our estimates would be reliable. That is, we would have to compute: 

m 

and is this quantity that could be expected to be zero. At this point two issues become 
apparent: 

i) Why is that so? why can this average be expected to yield zero but not each individual 
aim as in d?]). 

ii) Empirically, does this hold? In other words, is the actual average of observed complex 
quantities in ([8]), in fact, zero or is it not?. 



We will explore issues connected to these questions within the approach pioneered in 4| 
and which seems to have more potential to dealing with such questions than the standard 
one. 

The paper is organized as follows: In section II we review the standard picture for pri- 
mordial non-Gaussianities, in section III we review the collapse models description for the 
inflationary origin of the seeds of the cosmic structure. In section IV we focus on the statis- 
tical aspects of the primordial non-Gaussianities and propose new characterizations of the 
quantities associated with the bispectrum. Finally, in section V we discuss our findings. 
The conventions we will be using include a (— , +, +, +) signature for the space-time metric. 
We will use units where c = 1 but will keep the gravitational constant G and h explicit 
throughout the paper. 



II. THE STANDARD PICTURE FOR THE PRIMORDIAL NON-GAUSSIANITIES 

This section will briefly review the standard accounts on the primordial non-Gaussianities 



following closely references |14l-ll6j. There is absolutely no original work in this section; we 
simply present here the usual treatment on the subject following what is commonly found in 
the literature in order to compare with our own approach and discuss the main diferences. 
7or more details and derivations we refer the reader to the comprehensive review by Komatsu 

Ti n 

17|, Bartolo et al. 18| and the references cited therein. 

Historically, non-Gaussianity as a test of the accuracy of perturbation theory was first 
suggested by Allen et al. [19| . However must of its importance to date relies on the premise 
that it will play a leading role in furthering our understanding of two fundamental aspects 
of cosmology and astrophysics [20] : 

• The physics of the very early Universe that created the primordial seeds for large-scale 
structures, and 

• The subsequent growth of structures via gravitational instability and gas physics at 
later times. 

Within the standard approach, by non-Gaussianity, people refer to any small deviations 
in the observed fiuctuations from the random field of linear, Gaussian, curvature pertur- 
bations. The curvature perturbations, \l/, generate the CMB anisotropy, 5T /T. The linear 
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perturbation theory gives a linear relation between \l/ and 5T /T on large-scales (where the 
Sachs- Wolfe effect dominates) at the decoupling epoch, i.e. 5T /T ~ (1/3)^. It follows from 
the relation, (5T oc \E', that if \1/ is Gaussian, then 5T is Gaussian, but what exactly does one 
mean by Gaussian at the observational level?. 

One of the most important results of the inflationary paradigm is that the CMB 
anisotropy arises due to curvature perturbations which in turn are produced by quan- 
tum fluctuations. In the standard single field slow-roll scenario, these fluctuations are 
due to fluctuations of the inflaton fleld itself when it slowly rolls down its potential ^(0). 
Within this approach, the primordial perturbation is Gaussian; in other words its Fourier 
components are uncorrelated and have random phases. When inflation ends, the inflaton 
oscillates about the minimum of its potential and decays, thereby reheating the Universe. 

In the inflationary paradigm, the perturbations of the fleld 50 and the perturbations of 
the curvature \l/ are treated as standard quantum fleldCl evolving in a classical quasi-De 
Sitter background space-time. The quantity of observational interest is called the power 
spectrum of the curvature perturbation P^{k,ri). The power spectrum is obtained from 

(0|^(f,77)^(y,r7)|0), (9) 

where |0) is called the Bunch-Davies vacuum and represents the initial state of the fleld -0, 
which is the Mukhanov- Sasaki fleld variable deflned in (J2]). 

It is precisely at this step where a subtle issue arises, namely that in the standard picture 
one is given various and distinct arguments (e.g. decoherence, horizon crossing, many-worlds 
interpretation of QM, etc.) to accept the identiflcation 



(0|^(f,r/)^(y,r7)|0) = ^(f,r/)^(y,r7), (10) 

where \l/(x, r^) now stands as a classical stochastic fleld and the overline denotes the average 
over an ensemble of Universes. In other words, the value of the fleld \l/ in each point (x, rj) 
varies from each one of the members of the ensemble of "Universes", with a variance \&^. 
Therefore, the power spectrum Pm{k,rj) is deflned in terms of the Fourier components of 
^(f,?7) by 



In fact they are both part of a unified field v. 



^M^piv) ^ i2n)'5ik + k')P^ik,v). (11) 

Consequently, the power spectrum is related to the 2-point function through 

/fin SlTl nT" 

—V^{k,r])^—, (12) 

with r = \x — y\ and we also used the definition of the dimentionless power spectrum 
Vxi,{k,ri) = Pi!,{k,ri)k^/2iT'^. From the above expression one can find the variance ^^^ 



oo 



dk 



m\x,Tl) = j —V^{k,7]). (13) 

The expression (IT3l) diverges generically. In particular, we know that the spectrum of 
the primordial curvature perturbation is roughly P^,{k,vi) oc k~^. That is, V^iik^rj) is nearly 
constant (i.e. independent of fc); therefore (I13p diverges in a logarithm fashion for A; — )■ 



and A; — 7- oo. The way the standard pictures deals with this issue 21j is to establish a fcmax 
equal to the "horizon" , and work in a cubic box of physical size aL much larger than the 
Hubble radius. Thus, 



^ — = P^(r7)ln-^-^. (14) 

That is, in order to avoid the divergence in '^'^ one is forced to introduce some particular 
values of k as cut-offs. 

The question that arises now is how can we do an average over an ensemble of Universes 
if we have observational access to just one -our own- Universe. In the next subsection, we 
will show how the standard approach deals with this issue. In the following we will accept 
the validity of (TTOj). 

If \I'(x, rf) is Gaussiaio, then the two-point correlation function ([9]) specifies all the statis- 
tical properties of \E'(x, r^), for the two-point correlation function is the only parameter in a 
Gaussian distribution. If it is not Gaussian, then we need higher-order correlation functions 
to determine the statistical properties. 

For instance, a non- vanishing three-point functioql 

^ That is, there exist some physical mechanism for which the quantum variable ^{x, rj) becomes a classical 

stochastic field *I'(a;, 77) with Gaussian distribution. 
^ Similarly as the 2 point correlation function, the standard approach relies on the identification 



(0|*(f,77)*(y,77)*(z, 77)10) =*(f,77)*(y,?7)*(z,77). (15) 
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^{x,r])^{y,r])^{z,r]), (16) 

is an indicator of non-Gaussian features in the cosmological perturbations. The Fourier 
transform of the three-point function is called the bispectrum^ and is defined as 



^g^^,,^^,^^ = {27ry6{h + k2 + h)B^{ku k2, ks). (17) 

The importance of the bispectrum comes from the fact that it represents the lowest order 
statistics able to distinguish non-Gaussian from Gaussian perturbations. 

The delta function in (IT7j) enforces the triangle condition, that is, the constraint that the 
wavevectors in Fourier space must close to form a closed triangle, i.e. ki + k2 + k^ = 0. 
Different inflationary models predict maximal signal for different triangle configurations. 
The standard approach of the study of the structure of the bispectrum is usually done 
by plotting the magnitude of Bxi,{ki,k2,k^){k2lki)'^{k^/kiY (with \ki\ = ki) as a function 
of k2/ki and fcs/fci for a given fci, with a condition that ki > k2 > k^ is satisfied. The 
usual classification of various shapes of the triangles uses the following names: squeezed 
(ki ~ /c2 ^ ^3)5 elongated (fci = k2 + ^3), folded (fci = 2k2 = 2^^ ), isosceles {k2 = k^) and 



equilateral {ki = k2 = k^). Within the cosmology community |22N24| . these shapes of non- 
Gaussianity are potentially a powerful probe of the mechanism that creates the primordial 
perturbations. 

One of the first (and most popular) ways to parameterize non-Gaussianity phenomeno- 



logically was via a small non-linear correction to the linear Gaussian perturbation 25|, l26|] , 



^(f , 77) = ^l(^, v) + ^nl(^, v) 



Mx,v) + flWii^^v) - m^,v)], (18) 



where "^hix, 77) denotes a linear Gaussian part of the perturbation, and the variance "^f^ix, rj), 
is implemented in the same sense as presented in ( jT^ . Henceforth, we call f^l the lo- 
cal non-linear coupling parameter which determines the "strength" of the primordial non- 
Gaussianity. This parametrization of non-Gaussianity is local in real space and therefore is 



In the following we will not write the explicit dependance of the conformal time 77 unless it leads to possible 
confusion. 
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called local non-Gaussianity. In this local model the contributions from 'squeezed' triangles 
are dominant, that is, with e.g. fcs ^ ki,k2- Using f lTSj) and flT7|) the bispectrum of local 
non-Gaussianity may be derived 

In the standard picture, the non-Gaussianity produced by many single field slow-roll 
models, is considered small and likely unobservable. However, a large detectable amount of 



non-Gaussianity can be detected when any of the following conditions are violated 18 



27 



28|: 



• Single Field. There was only one quantum field responsible for driving inflation 

• Canonical Kinetic Energy. The kinetic energy of the quantum field is such that 
the speed of propagation of fluctuations is equal to the speed of light. 

• Slow Roll. The evolution of the field was always very slow compared to the Hubble 
time during inflation 

• Initial Vacuum State. The quantum field was in the preferred "Bunch-Davies vac- 
uum" state. 

It is expected that the primordial non-Gaussianity produced by inflation is very small, 
i.e. undetectable only when all of the above conditions are satisfied. 

A. Non-Gaussianity in the CMB 

In this subsection, we present the standard connection between the primordial bispectrum 
at the end of inflation and the observed bispectrum of CMB anisotropies. 

1. Theoretical predictions for the CMB bispectrum from inflation 

As we mentioned in section [U the temperature anisotropies are represented using the aim 
coefficients of a spherical harmonic decomposition of the celestial sphere, 
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^(^,<^) = ^a,„r,^(^,<^), (20) 

Im 

and the curvature perturbation \l/ is imprinted on the CMB multipoles aim by a convolution 
involving the called transfer functions A(/c) representing the linear perturbation evolution, 
through equation ([6]): 



The CMB bispectrum also called the angular bispectrum is defined as the three point 
correlator of the aim 



-"mim2m3 — Q^Zimi Q^i2m2 0^/3^3 • {^^) 

At this point, the standard picture lead us to another subtle issue, that is, the overline 
in fl^T]) denotes in principle an average over an ensemble of Universes. In reality, we cannot 
measure the ensemble average of the angular harmonic spectrum, but one realization such 
as {aiimi, ««2m2; • • • ;^m„}- To ovcrcome this issue, the standard approach relies on the 
ergodic assumption [21.]. The ergodicity of a system refers to that property of process by 
which the average value of a process characteristic measured over time is the same as the 
average value measured over the ensemble. If one accepts the common supposition that 
the inflationary perturbation is indeed ergodic, then one expects the volume average of the 
fluctuations to behave like the ensemble average: the Universe may contain regions where the 
fluctuation is atypical, but with high probability most regions contain fluctuations with root- 
mean-square amplitude close to a. Therefore the probability distribution on the ensemble, 
which is encoded in ( 12T]) . translate to a probability distribution on smoothed regions of a 
determined size within our own Universe. 

After the above analysis, we continue with the calculation relating the primordial bispec- 
trum with the angular bispectrum. By substituting Qj in (|2T|) . one obtains 
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- j / dkidk2dh {kik2hfB^{ki, k2, ks)A{ki)A{k2)A{k3) x 



ji^{kiRD)ji2ik2RD)ji3ihRD) / rfa; X ii^{kix)]i^{k2x)ii.^{k'ix) x 

/ dfia;>^imi(£)'i^i2m2(^)^3m3(^)5 (22) 

where in the last hne we have integrated over the angular parts of the three ki and used the 
exponential integral form for the delta function that appears in the bispectrum definition 
(fT7|) . The last integral over the angular part of x is known as the Gaunt integral, which can 
be expressed in terms of Wigner 3-j symbols as 



^/™"^' ^ j dn,Yi,„,,{x)Yi,^,{x)Y,, 



(2/i + l)(2/2 + l)(2/3 + l) ( h h h\ ( k h h 

mi 1712 ^^3 



47r 




(23) 



The fact that the bispectrum -8^1^^27113 consists of the Gaunt integral, GJ^lJl^"^^, implies 
that the bispectrum satisfies the triangle conditions and parity invariance: mi+m2+"^3 = 0, 
h + I2 + h = even, and \li — lj\ < Ik ^ k + Ij for all permutations of indices. 

One thus can write 

Tjhhh — nmim2m3i (0A\ 

^mim2mz ~ ^hhh "hhh^ \'^^) 

where ftziij^a is an arbitrary real symmetric function of li,l2 and l^. This form, fl2^ . is 
necessary and sufficient to construct generic B^^^^^^^ satisfying rotational invariance; thus 
in the literature one encounters h^i^i.^^ more frequently than B'^'-^^^^^. The quantity bi-^i^i^ is 
called the reduced bispectrum as it contains all the physical information in -Bmini2m3- Since 
the reduced bispectrum does not contain the Wigner 3-j symbol, which merely ensures the 
triangle conditions and parity invariance, it is easier to calculate physical properties of the 
theoretical bispectrum. 
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In the standard picture one assumes that if there is a non-trivial bispectrum then it has 
arisen through a physical process which is statistically isotropic^, so we can employ the 
angle- averaged bispectrum Bi-^i^i^ without loss of information, that is |l7l. Il8l|. 



Bhhh — 2^ I aUmi««2m2«Z3m3- (25) 

m, y mi m2 ms J 
We now can obtain a relation between the averaged bispectrum, Bi^^i^i^ and the reduced 
bispectrum bi-^i-^i^ by substituting f l24|) into ( 125|) . 



; (2/i + 1) (2/2 + 1) (2/3 + 1) hkh , 

BiM, = \l — 1 ' ''" 



where the identity. 



all m \ "^1 "^2 ^3 y V 47r \^ ' 

was used. The reduced bispectrum obtained from ( l22l) then takes the much simpler form 



^^1^2/3 = (-) / dhdk2dk3 {hk2h)^B^{ki,k2,k3)A{h)A{k2)A{h) x 

POD 

Jh{kiRD)ji2{k2RD)Jh{hRD) I dx x^ji,{kix)ji^{k2x)ji.,{k3x). (28) 

Jo 

This is the main equation for this section, since it explicitly relates the primordial bispec- 
trum, predicted by the standard inflationary theories, to the averaged bispectrum (through 
(^^) obtained from the CMB angular bispectrum ai^^^T^Ji^m^^^^sms- This formula is entirely 
analogous to the well known relation linking the primordial power spectrum Pq,{k) and the 
CMB angular power spectrum Q, i.e., 

Ci = - f k^P^ik)A^{k)jf{kRD)dk. (29) 



^ Although it would be interesting, and possibly a more realistic approach to the problem, to proceed in 
the analysis without this assumption. 
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2. Measuring 'primordial non-Gaussianity from the CMB 

As we mentioned before, in most inflationary models, the parameter characterizing pri- 
mordial non-Gaussianity is /nl- Thus the next task within the standard picture is to es- 
timate /nl from the CMB data-set. That is, one chooses the primordial model that one 
wants to test, characterizing it through its bispectrum shape. One then proceeds to esti- 
mate the corresponding amplitude f^^'^^^ from the data. If the final estimate is consistent 
with /^L^°' = 0, one concludes that no significant detection of the given shape is produced 
by the data, but one still determines important constraints on the allowed range of /nl'^*^^ 
Note that ideally one would like to do more than just constrain the overall amplitude, and 
reconstruct the entire shape from the data by measuring configurations of the bispectrum. 
However, the expected primordial signal is too small to allow the signal from a single bis- 
pectrum triangle to emerge over the noise. For this reason one studies the cumulative signal 
from all the configurations that are sensitive to /nl"^*^^- 

Given the above analysis, the standard picture then makes use of estimation theory to 
extract an estimate for /nl- An unbiased bispectrum-based minimum variance estimator for 



the nonlinearity parameter can be written as 29|, |30 | 



1 V-~\ / n ^2 '3 1 -D; ; ; 

/nl = ^^ \ 77r~Fr7r\ — ("/i»nia/2m2«i3m3)obs, (30) 

where -B^*] ^ is the angle averaged theoretical CMB bispectrum for the model in consideration 
with /j^L = 1) Q is the observed angular spectrum and aim are the multipoles of the 
observed CMB temperature fluctuations. The normalization A^ is calculated requiring that 
the estimator to be "unbiased", i.e. the averaged value is equal to the "true" value of 
the parameter, (/nl) = /nl- If the bispectrum Bi^i^i^ is calculated for /at^, = 1 then the 
normalization takes the following form 



, ChCi^Ci^ 

''i 

The estimator for non-Gaussianity (l30l) is then simplified using ( l28l) and (l26l) to yield 
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•^NL = ^X^ / dn£Yi^rni{x)Yi2^2{x)Yi^rn3{x) / X^cixj/, (fcix) jj^ (^23;) j/3 (/CsX) (C;^ ^C,^ ^C,^ ^)obs X 
hrrii 
/r,\ 3 ,. 



TT 



dkidk2dk3ikik2k3fB{h, k2, k3)A{ki)A{k2)A{k3)ji,{kiRD)jl,{k2RD)jls{k3RD) X 



where B{ki, k2, k^) is the primordial bispectrum obtained from the 3-point function as de- 
fined in ( 1T71) . In this manner the sought constraints are obtained. The best results, corre- 
sponding to the so called Jbcal, equilateral and orthogonal shape of non-Gaussianities using 



the WMAP 7-year data |3l| yield f^l""^ = 32 ± 21 (la), /^r' = 26 ± 140 (Ict) and f^^"^ 



= -202 ± 104 (la). 

III. THE COLLAPSE MODEL ACCOUNT FOR THE INFLATIONARY ORIGIN 
OF COSMIC STRUCTURE 

Before proceeding it seems worthwhile to briefly explain the view we take regarding 
quantum physics and Einstein's gravity. The framework we adopt is based on a description 
of the problem that allows, at the same time, the quantum treatment of other flelds and a 
classical treatment of gravitation. That is the realm of semi-classical gravity together with 
quantum fleld theory in curved space-time. We will assume to be a valid approximation 
for most of the time, with the exception associated precisely with the dynamical collapse 
as we will explain below. Such description of gravitation in interaction with quantum flelds 
is reflected in the semi-classical Einstein's equation: R^i, — {l/2)g^^R = 87rG{T^u), whereas 
the other flelds, including the inflaton, are treated in the standard quantum fleld theory 
fashion. It seems clear that this approximated description would break down in association 
with the quantum mechanical collapses or state jumps, that we consider to be part of the 
underlying quantum theory containing gravitation. The reason for this breakdown is simply 
that the left hand side of the previous equation does not contain any divergences, while the 
right hand side, and in particular its divergence 'Wp{T^^), will have discontinuities associated 
with the quantum jumps. 

In this setting we start from the assumption that, in accordance with the standard infla- 
tionary accounts, and as mentioned before, the state of the Universe before the time at which 
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the seeds of structure emerge is described by the H&I Bunch-Davies vacuum state for the 
matter degrees of freedom (D.O.F.) and the corresponding H&I classical Robertson- Walker 
space-time. 

Then, we assume that at a latter stage, the quantum state of the matter fields reaches 
a stage whereby the corresponding state for the gravitational D.O.F. is forbidden, and a 
quantum collapse of the matter field wave function is triggered by some unknown physical 
mechanism. In this manner the state resulting from the collapse of the quantum state of the 
matter fields needs not to share the symmetries of the initial state. After the collapse, the 
gravitational D.O.F. are assumed to be, once more, accurately described by Einstein's semi- 
classical equation. However as (Tfjn) for the new state, needs not have the symmetries of the 
pre-collapse state, we are led to a geometry that generically will no longer be homogeneous 
and isotropic. 

The starting point of the specific analysis is the same as the standard picture, i.e. the 
action of a scalar field coupled to gravity: 

S = j d'x,r9[-^R[g] - 1/2V.0V,0(?«^ - \/(0)], (33) 

where stands for the infiaton and V for the infiaton's potential. One then splits both, met- 
ric and scalar field into a spatially homogeneous ('background') part and an in-homogeneous 
part ('fluctuation'), i.e. g = Qq + 5g, (p = (pQ + 6(f). 

The background is taken to be the spatially flat Friedmann-Robertson Universe with 
line element ds"^ = a{riY [—drf + 6ijdx'^dx^], and the homogeneous scalar field 0o(^)- The 
evolution equations for this background are scalar field equations, 

• 2 

4>o + 2% + a%V[(P] = 0, 3^ = 47rG'(02 + 2aV[0o]). (34) 

(Jj (Jj 

The scale factor solution corresponding to the infiationary era of standard infiationary cos- 
mology, written using a conformal time, is: a{ri) = — l/[iJ|(l — e)ri] with iJ| ~ {8'k/3)GV, 
e = 1 — 'H/'H? is the slow-roll parameter which during the infiationary stage is considered 
to be very small e ^ 1, thus Hj ^ constant, and with the scalar 0o fi^ld in the slow roll 
regime, i.e. 0o = —{a?/3a)V'. According to the standard infiationary scenario, this era is 
followed by a reheating period in which the Universe is repopulated with ordinary matter 
fields, a regime that then evolves towards a standard hot big bang cosmology regime leading 
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up to the present cosmological time. The functional form of a{ri) during these latter peri- 
ods changes, but we will ignore those details because most of the change in the value of a 
occurs during the inflationary regime. We will set a = 1 at the "present cosmological time" , 
and assume that inflationary regime ends at a value oi rj = rjQ, negative and very small in 
absolute terms {rjo c^ —10^^^ Mpc). 

Next we consider the perturbations. We shall focus on the scalar perturbations and 
ignore for simplicity the gravitational waves. Working in the so called longitudinal gauge, 
the perturbed metric is written as: 

ds^ = a{r]f [-(1 + 2m)dif + (1 - 2m)5ijdx'dx^] , (35) 

where \E' stands for the scalar perturbation usually known as the Newtonian potential. 

The perturbation of the scalar field is related to a perturbation of the energy momentum 
tensor, and refiected into Einstein's equations which at lowest order lead to the following 
constraint equation for the Newtonian potential: 

V^^ = 47rG0o(50 = s50, (36) 

where we introduced the abbreviation s = 47iG(f)o. 

Now we consider in some detail the quantum theory of the field 6(f). It is convenient to 
work with the rescaled field variable y = a6(f) and its conjugate momentum tt = 6(f)/a. For 
simplicity we set the problem in a finite box of side L, which can be taken to oo at the end 
of all calculations. We decompose the field and momentum operators as 

k k 

— * 

where the sum is over the wave vectors k satisfying kiL = 27rnj for i = 1, 2, 3 with nj integer 
and where ^^(r/) = yk{^)d^ + yl{r])d ^ and %(?]) = gk{^)a^ + gl{r])a} -. with the usual choice 
of modes: 

Vkiv) = -7^ ( ^ ^ ~fc ) ^^P(-^^^)' 9k{v) = ^^Y 2 exp(-i/i;r7), (38) 



which leads to what is known as the Bunch-Davies vacuum. 

Note that according to the point of view we discussed at the beginning of this section, 
and having at this point the quantum theory for the relevant matter fields, the effects of 
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the quantum fields on the geometrical variables are codified in the semiclassical Einstein's 
equations. Thus equation ( 136|) must be replaced by 



V^^ = 47rG0o50 = s{6<j)) = is/a){n). (39) 

At this point one can clearly observe that if the state of the quantum field is in the vacuum 
state, the metric perturbations vanish and thus the space-time is homogeneous and isotropic. 

Our proposal is based on consideration of a self induced collapse which we take to operate 
in close analogy with a "measurement" (but evidently, with no external measuring apparatus 
or observer involved). This lead us to want to work with Hermitian operators, which in 
ordinary quantum mechanics are the ones susceptible of direct measurement. Therefore we 
must separate both ^^(r/) and vr^(r7) into their real and imaginary parts yj:{ri) = yjf^ijf) + 
Wk^iv) ^^^ ^fc(^) — ^k^iv) ^^^kiv) so that the operators y^' {rj) and tt^' {rj) are hermitian 
operators. 

So far we have proceeded in a manner similar to the standard one, except in that we 
are treating at the quantum level only the scalar field and not the metric fiuctuation. At 
this point it is worthwhile to emphasize that the vacuum state defined by a^^'''^|0) = is 
100% translational and rotationally invariant. That is, under spatial translations T{di) = 
exp[iPidi] and rotations Rx{Oi) = exp[iZ/(x)j^j], we have T{di)\0) = and Rx{Oi)\0) = 0, 
with Pi and L{x)i the linear and the angular momentum operators, and di and 6i parameters 
labeling the transformations. 

For the next step we must specify in more detail the modeling of the collapse. Then, 
take into account that after the collapse has taken place, we must consider the continuing 
evolution of the expectation values of the field variables until the end of infiation and even- 
tually up to the last scattering hypersurface (in fact, if we want to actually compare our 
analysis with the detailed observations, we must evolve also through the reheating period 
and through the decoupling era up to today's Universe. This, however, is normally taken 
into account through the use of appropriate transfer functions, and we will assume that the 
same procedure can be implemented after the present analysis). 

We will further assume that the collapse is somehow analogous to an imprecise measure- 
menlo of the operators jjj.' {r]) and n^' (77). Now we will specify the rules according to 



^° An imprecise measurement of an observable is one in which one does not end with an exact eigenstate 
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which collapse happens. Again, at this point our criteria will be simplicity and naturalness. 
What we have to describe is the state |6) after the collapse. It turns out that, for the goals 
at hand, all we need to specify is the quantity d^' = {Q\d^' |G), as this determines the 
expectation value of the field and momentum operator for the mode k at all times after the 
collapse. 

It is natural to assume that after the collapse, the expectation values of the field and 
momentum operators in each mode, will be related to the uncertainties of the pre-coUapse 
state (recall that the expectation values in the vacuum state are zero). In the vacuum state, 
yj: and tt^ individually are distributed according to Gaussian wave functions centered at 
with spread {Ay^)l and (A7r^)Q, respectively. 

We might consider various possibilities for the detailed form of this collapse. Thus, for 
their generic form, associated with the ideas above, we write: 



'k 



{yf\vl))e = Aix£\/(A^f )g = Xix^';Mvl)\^/hLy2, (40) 



{-.'" 



\vt))e = A2x£y(Avrf)2, = \,x^,;i\g,{r^l)\y^h!Jj2, (41) 

where x^. \ , x^ 2 are selected randomly from within a Gaussian distribution centered at zero 
with spread one and rj'^ represents the time of collapse for each mode. Here Ai and A2 are 
two real numbers (usually or 1) that allow us to specify the collapse proposal we want 
to consider. At this point, we must emphasize that our Universe corresponds to a single 
realization of these random variables, and thus each of these quantities x^ \ , x^ 2 has a single 
specific value. On the other hand, we will be using the values for Ai and A2 to characterize 
the different collapse schemes, e.g.: i) Ai = 0, A2 = 1, ii) Ai = A2 = 1 (which we call the 
"symmetric scheme"). It is clear that one can devise many other models of collapse, a good 
fraction of them can be described within in the scheme above, while others require a slightly 
modified treatment [7^. Still, there are surely many other possibilities which we have not 
even thought about and which might require drastically modified formalisms. 

Finally for each model we obtain the information giving the relevant expectation values of 
the field operators in the post collapse state |G). That is, from the equations above, one can 
of that observable but rather with a state which is only peaked around the eigenvalue. Thus, we could 
consider measuring a certain particle's position and momentuni so as to end up with a state that is a 
wave packet with both position and momentum defined to a limited extent, and which, of course, does 
not entail a conflict with Heisenberg's uncertaint^-ibound. 



solve in each case for the quantities d-' , and using the result in the evolution equations for 
the expectation values (i.e. using Ehrenfest's Theorem) one obtains {y^^'^irj)) and (vr^'^'^(^)) 
for the state that resulted from the collapse, for all later times. The explicit expressions for 
the {yf\ri))Q and (7r5^(r/))0 are 



{yf\'n))e 



COsDfc / 1 

kr] 



k 



Zk 



k 



kr]Zk 



+ i™sfl.-5i^)(y?'(,a) 



G, 



{-ri^^)e 



(42) 



(^?''(^))e 



cos Dfc + 



sin Di 



Zk 



where D^ = ki] — Zk and z^ = kr]% 



■ R,I 



RJ I 



n^'\vd)e-ksmD,{y^'\7^l))e, 



(43) 



With this information at 
the collapse of all the mode; 



land we can now compute the perturbations of the metric after 



fl 



A. Connection to Observations 

Now, we must put together our semi-classical description of the gravitational D.O.F. and 
the quantum mechanics description of the infiaton field. We recall that this entails the semi- 
classical version of the perturbed Einstein's equation that, in our case, leads to equation 
(j39|) . The Fourier components at the conformal time rj are given by: 



*'<"' = -i/Ii^**''""' 



(44) 



where we used that during inflation s = AnGcpo = \Je/2aHi/Mp^ with Mp the reduced 
Planck's mass Mp = h^/^SnG). The expectation value depends on the state of the quantum 
field, therefore, as we already noted, prior to the collapse, we have ^E'^(?7) = 0, and the 
space-time is still homogeneous and isotropic at the corresponding scale. However after 
the collapse takes place, the state of the field is a different state with new expectation 
values which generically will not vanish, indicating that after this time the Universe becomes 



^^ In fact, we need only be concerned with the relevant modes, those that affect the observational quantities 
in a relevant way. Modes that have wavelengths that are either too large or too small are irrelevant in 
this sense. 
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anisotropic and in- homogeneous at the corresponding scale. We now can reconstruct the 
space-time value of the Newtonian potential using 

^(^'^) = X5E e^'%^(^)' (45) 

k 

to extract the quantities of observational interest. 

In order to make contact with the observations we shall relate the expression (j33]) for 
the evolution of the Newtonian potential during the early phase of accelerated expansion, 
to the small anisotropics observed in the temperature of the cosmic microwave background 
radiation, 6T{9,ip)/To with Tq ^ 2.725i^ the temperature average. They are considered as 
the fingerprints of the small perturbations pervading the Universe at the time of decoupling, 
and undoubtedly any model for the origin of the seeds of cosmic structure should account 
for them. As already mentioned in section [H these data can be described in terms of the 
coefficients aim of the multipolar series expansion, i.e. equation (E]). The different multipole 
numbers / correspond to different angular scales; low / to large scales and high / to small 
scales. At large angular scales (/ < 20) the anisotropics in the CMB arise due to the Sachs- 
Wolfe effect. That effect relates the anisotropics in the temperature observed today on the 
celestial sphere to the inhomogeneities in the Newtonian potential on the last scattering 
surface, 

^Z^e,ip) = \^iVD,XD). (46) 

Here rjD is the conformal time of decoupling which lies in the matter-dominated epoch, 
and xd = -Rz)(sin 9 sin ip, sin 9 cos ip, cos 9), with Rd the radius of the last scattering surface. 
Furthermore, using fllS]) and e*^'^-° = 47r ^;^i'ji(A;_R/))Yim(6', v5)y;^(A;), the expression (^ 
for aim can be rewritten in the form i^. The transfer function A(A;) represents the evolution 
of the Newtonian potential from the end of inflation rjn to the last scattering surface at the 
time of decoupling rjD, i.e. "^j^^tid) = A{k)'^j:{riji). 
Substituting (gS]) in (gl]) and using (gO]), (gl]) gives 



*,(,«) ^ -<^''"'^^' 



2V2MpA;3/2 
Finally using, (g7|) in ([6]) yields 



Aal cosDfcH — ^ j(xf2+«a;fc2)--^isinL'fc( 1+— j (a^Ji+^a^li^ 



(47) 
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[^k,2 + ^^A:,2J 



"- = 'WW^ ^ A(A;),K^/?,)1^:,(^) [a, ( cos D, + ^) , 

(48) 



k 

- AisinDfcf 1 + — j ixj}i +ixi^^] 

note that in ( HTll and ( HHll . D^ is evaluated at r^ij, i.e. Dk{r]R) = krjR — Zk- 

It is worthwhile to mention that the relation for aim with the Newtonian potential, as 

obtained in (HHl) within the collapse framework has no analogue in the usual treatments of 

the subject. It provides us with a clear identification of the aspects of the analysis where the 

"randomness" is located. In this case, it resides in the randomly selected values x^^ \, Xf, \ that 

appear in the expressions of the collapses associated with of each modes. Here, we also find 

a clarification of how, in spite of the intrinsic randomness, can we make any prediction at all. 

The individual complex quantities aim correspond to large sums of complex contributions, 

each one having a certain randomness but leading in combination to a characteristic value in 

just the same way as a random walk made of multiple steps. In other words, the justification 

for the use of statistics in our approach is: The quantity aim is the sum of contributions 

from the collection of modes, each contribution being a random number leading to what is 

in effect a sort of "two-dimensional random walk" whose total displacement corresponds to 

the observational quantity. Nothing like this can be found in the most popular accounts, 

in which the issues we have been focussing on, are hidden in a maze of often unspecified 

n 

assumptions and unjustified identifications [5,]. 



Thus, according to f HHj) all the modes contribute to aim-, with a complex number. If 
we had the outcomes characterizing each of the individual collapses we could of course 
predict the value of this quantities. However, we have at this point no other access to such 
information than the observational quantities aim themselves. 

We hope to be able to say something about these but doing so requires the consideration 
of further hypothesis regarding the statistical aspects of the physics behind the collapse as 
well as the conditions previous to them. 

As is in general the case with random walks, one can not hope to estimate the direction of 
the final displacement, however one might say something about its estimated magnitude. It 
is for that reason we will be focussing on estimating the most likely value of the magnitude: 
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Win.? = ^jjJ^A{k)A{k%{kRn)Uk'Rn)Y;^{k)YiUk')^,^{VR)'^g*{VR)- (49) 

k,k' 

We can progress for instance by making some assumption allowing us to regard the specific 
outcomes characterizing our Universe as a typical member of some hypothetical ensemble of 
Universes. 

For example, we are interested in estimating the most likely value of the magnitude of 
|a;mP above, and in such hypothetical ensemble we might hope that it comes very close 
to our single sample. It is worthwhile emphasizing that for each / and m we have one 
single complex number characterizing the actual observations (and thus the real Universe 
we inhabit). For a given / for instance, we should avoid confusing ensemble averages with 
averages of such quantities over the 2/ + 1 values of m. The ensemble are figments of our 
imagination and there is nothing in our theories that would indicate that they are real. 



We can simplify things even further by taking the ensemble average ja^mp (the bar will 
from now on means that we are taking the ensemble average) and identifying it with the 
most likely value of the quantity, and needless is to say that these two notions are different 
in many types of ensembles. However, let us for the moment ignore this issue and assume 
the identity of those two values and look at the ensemble average of the quantity |«;mP 
which is given by 

k,k' 

One can for instance assume that the collapsing events are all uncorrelated, and then 
consider estimating the most likely value thus 



2 167r2 



c^imltAL = ^jj}_^A{k)A{k%{kRn)Mk'Rn)Y;^{k)YiUk')^dVR)'^k*iVR)- (51) 



k,k' 



Under the assumption of the validity of such approximation and the additional assump- 
tion that the random variables xj^^, xl i, xf^g^ ^i 2 ^^^ all uncorrelated, we obtain that all the 
information regarding the "self-collapsing" model will be codified in the quantity: 



^fc(%)^p(^«)- (52) 
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Furthermore, we can take the hniit —krjn — )■ in (152|) . which can be expected to be ap- 
propriate when restricting interested on the modes that are "outside the horizon" at the end 
of inflation, since these are the modes that give a major contribution to the observationally 
relevant quantities. 

Under those conditions, and with the help of (H7|) . after taking the continuum limit 

(L — i- oo), the quantity of interest becomes: 

hhH] f dk 



where some of 



function C{k 



£ 



Wim?M^. = 1^ y ^A\k)jf{kRn)C{k), (53) 

;he information regarding that a collapse has occurred is contained in the 



. The explicit form of C{k) is 

C{k) = Xlfl + L\ sin2 zk + \l( cos Zk - ^] . (54) 

As we have noted in previous works, this quantity becomes a simple constant if the collapse 
time happens to follow a particular pattern in which the time of collapse of the mode k 
is given by rjl = Z/k with Z a constant. In fact, the standard answer would correspond 
to C{k) = constant (which can be thought as an equivalent "nearly scale invariant power 
spectrum"). Thus, the result obtained for the relation between the time of collapse and the 
mode's frequency, i.e. ?7|A;= constant is a rather strong conclusion which could represent 
relevant information about whatever the mechanism of collapse is. A preliminary study of 
the effects of small deviations from such pattern for the "symmetric scheme" (Ai = A2 = 1); 
for the "Newtonian scheme" (Ai = 0, A2 = 1) and for a third scheme that falls outside the 
category contemplated here, have been carried out in [7|. 

It is quite clear that if the time of collapse of each mode do not adjust exactly to the 

pattern krjl = Z, then the various schemes of collapse, characterized by the values of Ai, A2 

or some other function C{k), would lead to different predictions for the exact form of the 

^^ The standard amplitude for the spectrum is usually presented as ex V/{eMp) oc H'j /{Mpt). The fact 

that e is in the denominator leads, in the standard picture, to a constraint scale for V . However, in 

(|53p e is in the numerator. This is because we have not used (and in fact we will not) explicitly the 

transfer function A(fc). In the standard literature it is common to find the power spectrum for the 

quantity Q{x), a field representing the curvature perturbation in the co-moving gauge. This quantity 

is constant for modes "outside the horizon" (irrespectively of the cosmological epoch), thus it avoids 

the use of the transfer function. The quantity C, can be defined in terms of the Newtonian potential as 

C = * + (2/3)('H"i* + *)/(l + w), with cj = p/p. For large-scale modes C ^ ^[(2/3)(l + uj)-^ + 1], and 
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during infiation 1 + a; = (2/3)e. For these modes 4 ~ '^ /e and the power spectrum is Vc^{k) — V^{k)/e^ ex 
H?/(M?^f) nc V/(fM'k). which contains the correct a^mnlitude. For a, detailed discussion reffaTdinff the 



spectrum, and comparing these predictions with the observations can help us to discriminate 
between the distinct collapse schemes. 

We end this section by noting that the treatment of the statistical aspects in the collapse 
proposal is quite different from the standard inflationary paradigm. We will deepen this 
discussion in the next section. However, at this point the differences should be evident. 
In the standard accounts, one is going from quantum correlation functions to classical n- 
point functions averaged over an ensemble of Universes; then one goes to n-point correlation 
functions averaged over different regions of our own Universe, and finally one relates this 
last quantity with the observable |a;mP- These series of steps are not at all direct and they 
involve a lot of subtle issues that the standard picture does not provide in a transparent way. 
On the other hand, in the collapse proposal, the observable la^mP is related to the random 
variables x^s through a "random walk", as we mentioned, the value of |a«mP corresponds to 
the "length" of the random walk. This random walk is associated to a particular realization 
of a physical quantum process (i.e. the collapse of the inflaton's wave-function), since we 
have only access to one realization -the random walk corresponding to our own Universe- 
the most natural assumption (but certainly not the only one) is that the average value of 



the length of the possible random walks, which corresponds to la^mP, is equal to the most 
likely value, i.e. to |tt«m|ML) ^"^^ ^^is in turn, is associated with the |Q;;mP of our observable 
Universe. 

IV. STATISTICAL ASPECTS 

The first thing we should now note is that there are several statistical issues at play, and 
that within our approach, various novel ones emerge. One aspect is the exact nature of the 
state previous to all collapses, i.e. the state characterizing the first stages of the inflationary 
regime and normally taken to be the Bunch-Davies vacuum. There are various aspects that 
might affect and modify the nature of that state: For instance, if the field is not truly a 
free field and self interactions are important one might find correlations between the various 
modes of the field. These effects could be manifest, for instance by non vanishing values of 
quantities like (0|^g^g,|0). However, we should be aware, not only of the inherent problems 
of accessing these associated with the fact that we have at our disposal a single Universe, 
but also that our Universe, including the relevant perturbations is not characterized by the 
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vacuum state, but rather by the state that results after the collapses of all the modes, and, 
it is quite clear that the collapse process itself can be a source of unexpected correlations. 
These would manifest themselves, for example, in correlations between the values taken by 
the xg's appearing in the collapse process and which we have so far assumed were different 
and independent quantities for each mode. 

Moreover, we have to note that the quantities that are more or less directly accessible 
to observational investigation are not the (G|^^|0) and the n-point functions in general, for 
the post-collapse state, but the various a/m's, and the latter are related to the former, as 
can be seen in f HHj) in a nontrivial way. In fact as we saw each aim corresponds to a sort 
of two dimensional random walk (i.e. a sum of complex quantities) and each of the steps 
is related to (6|y^|0). It is thus clear that there might be correlations between the various 
aimS simply due to the fact that they arise from different combinations of the same random 
variables. Of course, we should note that the particular analysis behind our analysis was 
based on the assumption that the elementary process was associated with the collapse of the 
observables y^ and their conjugate momenta according to (HOj) . (HT]) . It is clearly conceivable 
that the elementary process might have been associated instead with other observables. One 
simple possibility for those alternative observables, are the various options offered by linear 
combinations of the former. 

A. The new outlook on non-Gaussianities 

In this section we discuss the aspects that need modification in the study of primordial 
non-Gaussianities in view of the approach we have been discussing to the origin of the 
primordial fluctuations. 

The first point we should stress is that from the two aspects of cosmology mentioned in 
section HI we have seen that we have had to modified the first, namely the nature of the 
quantum state in order to be compatible with the existence, at the fundamental (quantum) 
level, of the inhomogeneities and anisotropics that are behind the emergence of structure 
and thus of everything -including observers- in our Universe. 

Namely, the standard physics of the very early Universe had to be supplemented with 
the collapse hypothesis in order to fully account for the process that created the primordial 
seeds for large-scale structure. Otherwise, we could not really identify the process by which 
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the inhomogeneity and anisotropies emerged from the initial vacuum. 

As in the standard approach, we take the curvature perturbations \I/ to be the generators 
of the CMB anisotropy, 5T /T, however, in our approach, the observed fluctuations are 
determined not just by the initial vacuum state, which is and remains homogeneous and 
isotropic, but also by the characteristics of the collapse process, besides, of course the effects 
of the late time physics. 

In this more precise and detailed approach, it is clear that even if the primordial state 
can be consider as Gaussian; in the sense that the corresponding n-point functions are 
completely determined by the 2-point functions, -and thus the odd n, ra-point functions 
vanish- it might still be possible for collapse process to drastically affect and modify this. In 
other words there exist in principle the possibility that the collapse process itself introduces 
non-Gaussian characteristics into the state. We will not discuss this possibility here but 
only point it out as something to have in mind, and as a topic for future research. 

As we have argued, the quantity of observational interest is not really (0|^(x, ?7)\l/(y, ?7)|0) 
as the argument to justify that in the standard approach, depends not only on accepting the 



identification (0|\E'(x, r7)\E'(y, ?7)|0) = '^(x,ri)'^{y,ri), where '^{x,ri) is taken to be a classical 
stochastic field and the overline denoting the average over an ensemble of Universes, but 
also, on a series of arguments indicating one can replace the ensemble averages with suitable 
spatial averages of quantities in our Universe. 

As a matter of fact, a clear example of how a careless approach to the statistics at hand 
can lead to wrong conclusions, is brought by the variance \E'^. We mentioned in section [III 
that \E'^ diverges generically if we do not introduce an ad-hoc cut-off for k. Therefore, if we 
consider the temperature fluctuations in a particular point xq of the CMB 2-sphere and we 
estimate it in terms of (0|\l/^(xo,''7)|0), we obtain a divergent quantity, but clearly, from the 
observational data, we know that these fluctuations of the mean temperature in a particular 
point are rather small ~ 10~^. Meanwhile, in the collapse proposal, we do not need to deal 
with these issues because the variables subjected to the collapse are not y{x, rj), n{x, rf) but 
the field modes yjiirj), T^^irf), i.e. the collapse does not occur in the position space, instead we 
consider it as a collapse for each mode k. Therefore, the quantities of observational interest, 
namely the la^mP's depend on the expectation values {y^ij]))^^ (%(^))0 i^i the state |6) 
after the collapse. 

As we saw in the Introduction, if we really took \E'(x, rj) to be Gaussian and allowed the 
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identification of its n-point functions with the observations, we would have to accept that 
such identification holds in particular for the 1-point function, and that would lead us to a 
clear conflict between theory and observation. 

Similarly we must be careful when we consider the three-point function: 



^(x,r/)^(y,r7)^(z,r/), (55) 

and use the identification with ensemble averages with the measured quantities as an indi- 
cator of non-Gaussian features in the cosmological perturbations. 



As we saw the bispectrum \E'^ vp^ vp^ = (27r)^(5(A;i + k2 + k^)Bxi,{ki^k2,kz) is usually 
said to represent the lowest order statistics able to distinguish non-Gaussian from Gaussian 
perturbations, because Gaussianity is identified with the requirement that all statistical 
information is contained in the 2-point functions, and thus implicitly with the vanishing of 
all ra-point functions with n odd. However, the lowest odd integer is 1 not 3, and as we 
already saw the serious issue that arises with the 1-point function, forces us to question the 
standard arguments. 

In fact looking anew at the quantities normally associated with the 1-point function we 
see that we have at our disposal not only the average quantities Ci, but also, for every value 
of / and m, the individual quantities aim- Each one of those correspond, in our approach, 
to different random walks. It could prove very interesting to study the distribution of these 
numbers themselves: namely we can look at the plot of, say the magnitudes \aim\ for a given 
value of /. This set of 2/ + 1 numbers can naturally be expected to display a Gaussian shape. 

This seems to be a particularly relevant analysis and it is not clear to us if something 
like that has been studied in the literature. It seems to us, that the traditional approach 
would have naturally lead to the consideration of that issue. Looking at the distribution 
of the corresponding phases should be equally enlightening. Moreover, as we mentioned in 
the discussion around the equation ([H]), it would be interesting to evaluate the quantity ai 
defined there, and compare the result with the any of the natural estimates for its value, 
particularly the expected ensemble average of its magnitude. 

Another point worth to note is, as we have mentioned already, it is usually believed 
that a large detectable amount of non- Gaussianity is to be expected when the initial state 
of the quantum field is not the preferred Bunch-Davies vacuum state. Nevertheless, in 
the collapse proposal we are considering, the quantum state of the field after the collapse 
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is 1 6) 7^ |0) (th e analysis of a particular characterization of the post-colapse state has 
been done in 32j). Therefore, the curvature perturbation responsible for the temperature 
anisotropics in the CMB, is due to the expectation values (^^)e and (7r^)0 which in principle 
could generate detectable non-Gaussianities, as these quantities are never considered in the 
standard accounts. Of course, a further exploration of these ideas is required. 

The other delicate issue related to the statistical aspects of the traditional approach is 
encoded in the ergodicity assumption which we will discuss in the following. For instance, 
as we already saw in section [TTl the CMB bispectrum was defined as the 3-point correlator 
of the aim through -Bmini2m3 — f^i^mTf^z^m^f^m^- The standard picture forces one to deal 
with the issue that the RHS represents an average over an ensemble of Universes, while 
we have but one realization {a^^mi, Cii2m2, ■ ■ ■ , Cii„m„}- To overcome this issue, the standard 
approach relies on an ergodicity assumption, and which identifies the average value of a 
process characteristic measured over time to be the same as the average value measured 
over the ensemble. 

There are various issues that lead one to be concerned about this assumption and the 
application to the situation at hand: The first thing we must be aware is that ergodicity 
is a property of systems in equilibrium and it is rather unclear why this should be valid 
regarding the conditions associated with the inflationary regime. 

Next, as we mentioned, the ergodicity assumption is translated in the case at hand into 
the notion that the volume average of the fluctuations behave like the ensemble average 
"the Universe may contain regions where the fluctuation is atypical, but with high probability 
most regions contain fluctuations with root-mean-square amplitude close to a " , and thus one 
argues that the probability distribution on the ensemble, translates to a probability distribution 
on smoothed regions of a determined size within our own Universe. 

There are three issues here: 

i) How do we go from the arguments supporting ergodicity in time averages, to corre- 
sponding arguments for spatial averages?. 

ii) Regarding the CMB, we in fact do not have access to the spatial sections that would 
allow us to investigate the space averages. We only have access to the particular 
intersection of our past light come with the 3-d hypersurface of last scattering: that is 
to a 2-sphere that we see as the source of the CMB photons that reach us today. Lets 
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call it from now on the CMB 2-spliere. \}_ 

iii) The third problem is that each one of the quantities of interest aim is itself already 
a weighted average over the CMB 2-sphere (with the weight function given by the 
corresponding Yim{0,f))- Therefore, what would be this new average we would be 
talking about in the arguments above?. In other words, if one is willing to accept that 
the ensemble averages should coincide with averages over the 2-sphere, why would one 
not also accept that the weighted average over the two sphere, should coincide with 
the equally weighted average over ensembles?. The latter would be zero so this is 
clearly unacceptable. 

It seems clear that we are dealing here with an orientation average: The different aim 
would mix among themselves if we were to redefine the orientation of the coordinate chart 
used to describe the CMB 2-sphere. When we look at the averages that are actually per- 
formed in connection with the study of the primordial spectrum we see these are indeed 
orientation averages. For instance, the observational quantity Cf^^ = ^^^^ J2m \(^im\'^ is just 
the orientation average value of the magnitude of the a^m's for a fixed value of /. In the same 
fashion, we see that the angle- averaged bispectrum -B^i./j./a (ES]), is an orientation average 
for fixed Vs and, as for the same reason as the 1-point function, it is quite unclear how to 
identify orientation averages with ensemble averages. Thus, the statistical analysis would 
be more transparent if one would focus on the distribution of the quantities B^^^^^^^^. 

As we saw, it is customary to take as an estimator for the nonlinearity parameter the 
quantity /nl defined in f l32|) . This seems a bit problematic as it involves a mixture of 
theoretical and observational quantities. Ideally one would like to have the two aspects 
rather well separated. In fact even within the standard approach, for the case of the 2-point 
functions we have on the one hand: the theoretical quantity 

Cf^ = - f k^P^{k)A^{k)jf{kRD)dk, (56) 

and on the other hand the observational quantity. 



^•^ We note in relation to this point that there are intrinsic problems in considering ergodicity of processes 

33I. 



within a two sphere as discussed in 
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This independence of the definitions allows one to cleanly compare theory and obser- 
vation. It thus seems that one would want to consider studying the aspects tied to non- 
Gaussianity, using a quantity that can be equally susceptible to theoretical and observational 
determination. Here we would like to propose, based on the considerations we have been 
discussing, the option we discuss below. 

First, motivated by the quantity defined in ( 125|) let us introduce the definition of the 
observed bispectrum as the orientation average 

^hllh — Z-^\ {Oihmiai2m2<^l3m3)ohs, (58) 

^^ y nil rn2 m^ J 
and the definition of the normalized observational reduced bispectrum as the quantity: 



K 



pbs 



hh 



-\ -1 



(2/1 + 1) (2/2 + 1) (2/3 + 1) / k h h 



SfZis^ (59) 



47^ I 

and finally let us define the magnitude of the bispectral fiuctuations as: 

•^«°ii2i3 — (2/ I l)(2l I l)(2l I 1) ^ M'^hmiOihmiO^hmsT ' " ^hhTa^'^^Khhi ' (^0) 

TTli 

One can then compare this pure observational quantity with the corresponding theoretical 
estimation corresponding to suitable ensemble average where each element of the ensemble is 
characterized by the specific random choice of the numbers xg that characterize the collapses. 

This seems to offer an approach to the issue at hand that indeed has the advantage 
of allowing a direct comparison between the purely observational quantities untainted by 
theoretical models, and the quantities that are purely defined in terms of such theoretical 
analysis. This in fact seems to share some of the spirit of the analysis made in references 
|l4J . Il5| . although our proposal provides a clear potion to compute the observational and 
theoretical quantities in complete separation, and that seems not to be available in the 
former. The reason for this seems easy to understand: The fact that we maintain a clear 
distinction between ensemble averages and orientation averages avoids the possibility of the 
confusion associated with the fact that the ensemble average of the quantity 
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appearing in (160!) vanishes identically. 

The detailed analysis of estimators like this will be carried out in future works, but we 
wanted to present it as an example of the type of analysis that is motivated by our approach 
to the whole question of emergence of structure from quantum fluctuations in the inflationary 
early Universe. 

V. PREDICTIONS AND DISCUSSION 

Focussing in trying to understand the essence of the emergence of inhomogeneous and 
anisotropic features from a quantum state that is homogeneous and isotropic and in the 
absence of a measurement procesO, has lead us to consider modifying the standard approach 
through the incorporation of the collapse hypothesis. 

We have seen in previous works that despite the fact that the motivation for such consid- 
erations seems to be purely philosophical and tied to issues like the measurement problem in 
quantum mechanics, the analysis have led to expect certain departures that could potentially 
be of observational significance. 

In previous works, we have focussed on two main observationally related issues: The 
shape of the spectrum and the question of tensor modes. We have argued previously that 
it would be very unlikely that one could find a scheme in which the function C{k) would be 
exactly a constant, and that some dependence of k will remain in any reasonable collapse 
scheme, simply because we do not expect those collapses to follow exactly the r]^ = Z/k rule 
for the time of collapse for each mode. These dependences lead to slight deviations from 
the standard form of the spectrum. In fact a preliminary analysis of this issue has been 
carried |7| , confirming these expectations. A more detailed analysis of the exact form of the 
spectrum incorporating the late time physics is needed, and under way, to study the traces 
left by such effects in the observational CMB spectrum, and as we have indicated, these can 
be searched for observationally. 

We have also stated, in previous works, that the most clear prediction of the novel 

paradigm we have been proposing, is the absence of tensor modes, or at least their very strong 

suppression. The reason for this can be understood by considering the semiclassical version 

^^ In the early Universe there are no observers or measuring devices, and in fact, the conditions for their 

emergence is the result from the breakdown of such symmetries, so it would seem very odd if one took 

a view that they are part of the cause of that breakdown. 
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of Einstein's equations and its role in describing the manner in which the inhomogeneities 
and anisotropies in the metric arise. As we have explained in our approach, the metric 
is taken to be an effective description of the gravitational D.O.F., in the classical regime, 
and not as the fundamental D.O.F. susceptible to be described at the quantum level. It 
is thus the matter degrees of freedom (which in the present context are represented by the 
infiaton field) the ones that are described quantum mechanically and which, as a result of a 
fundamental aspect of gravitation at the quantum level, undergo effective quantum collapse 
(the reader should recall that our point of view is that gravitation at the quantum level 
will be drastically different from standard quantum theories, and that, in particular, it will 
not involve universal unitary evolution). This collapse of the quantum state of the infiaton 
field leads to a nontrivial value for (T^u), which then generates the metric fluctuations. 
The point is that the energy momentum tensor contains linear and quadratic terms in the 
expectation values of the quantum matter field fluctuations, which are the source terms 
determining the geometric perturbations. In the case of the scalar perturbations, we have 
first order contributions proportional to 0o(<^0) while no similar first order terms appear as 
source of the tensor perturbations (i.e. of the gravitational waves). Of course, it is possible 



that the collapse scheme works at the leve 
metric fluctuations as has been presented in 



of the simultaneously quantized matter and 



34( 1 although, as explained there, we would find 



much harder to reconcile that with the broad general picture that underlies of our current 
understanding of physics. 

In the present work we have focused on the modified statistical considerations associated 
with this novel paradigm. We have argued that the collapse process itself could be the 
source of non-Gaussian features. We discussed some difficulties associated with the usual 
identification of measuring quantities with the quantum ra-point functions and particularly 
found that extending the standard arguments to the 1-point functions lead to disastrous 
disagreements with observations. 

We have shown that our approaches provides expressions which have no parallel in the 
standard formulations and which allow a precise identification of the location of the random- 
ness, as exemplified by our theoretical formula ( HHj) for aim in terms of the random numbers 

D r 

characterizing the collapses, namely the quantities x^' . This kind of expression facilitates 
all resulting statistical considerations and in particular, it is the basis for the theoretical 
estimation of the quantity ( 160|) . 
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We have proposed various novel ways to look into the statistical aspects of the problem: 

i) We indicate the importance of exploring the true nature of the one point function 
by studying the degree of deviation from zero from the complex quantity a°^^ = 

21+1 Z^Z ^Im ■ 



ii) We have argued that it is worthwhile to study the specific form of the distribution of 
the values of the observed quantities \a°^\ for each fixed /. 

iii) We have proposed new characterizations of the quantities normally associated with 
the bispectrum and the quantum 3-point functions which can be computed both, in 
purely theoretically, and in a completely observational fashion. This is the quantity 
defined in (l60l) . 

It is clear that this work represents only the first step in the study of the statistical aspects 
of the cosmic structure and its generating process during inflation, within the context of the 
new paradigm which centers on the collapse hypothesis. Much more work remains to be 
done, but we hope this can become a research avenue of great richness, and one which 
would lead to important insights, with possible implications not only for the generation 
of structure itself but for the modification of quantum theory, which would underly the 
collapse mechanism and which, as has been argued before, might have deeper origins at the 



quantum/gravity interface p, llO|, ll2, l35[ |. 
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